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ABSTRACT Accreting black holes and neutron stars in their hard (low) state show not only 
very similar X/'y-ray spectra but also that the behaviour of their light curves is quite similar 
which can be quantified as having similar power-density spectra and Fourier-frequency-dependent 
time/phase lags. Taken together this argues for a common mechanism of the X/7-ray production 
in these objects. This mechanism is probably a property of the accretion flow only since it does 
not depend on the nature of the compact object. In this paper, I review the observational data 
paying most attention to the properties of the temporal variability such as the time/phase lags 
that hopefully can help us to discriminate between different theoretical models. I also discuss 
the models developed to account for the basic observational facts. Particularly, I show that the 
commonly used Compton cloud models with constant temperature cannot explain variable sources 
without violating the energy conservation law. Alternative models where time lags are related to 
the spectral evolution during X-ray flares are discussed and compared with observations. Compton 
reflection from the outer edge of the accretion disc is shown to markedly affect the time lag Fourier 
spectrum. 

KEYWORDS: accretion, accretion discs; black hole physics; stars: neutron; stars: flare; stars: 
individual (Cygnus X-l); X-ray: stars. 

1. MODELS FOR THE FORMATION OF X/7-RAY SPECTRA 

X-ray and gamma-ray spectra of accreting black holes and neutron stars are de- 
convolved into (at least) two components: a soft component interpreted as emission 
from an optically thick accretion disc, and a hard tail associated with a hot (10-100 
keV) "corona" . Reviews of the spectral properties of Galactic black hole candidates 
(GBHs) can be found in Gilfanov et al. (1995), Tanaka & Lewin (1995), Grebenev 
et al. (1993, 1997), Grove et al. (1997), and Poutanen (1998). X/7-ray properties 
of radio-quiet active galactic nuclei (AGN) are reviewed by Zdziarski et al. (1997), 
Johnson et al. (1997), and Zdziarski (1999). Recent results on the broad-band spec- 
tra of accreting neutron stars are presented by Barret et al. (2000). 

An amusing fact is that super-massive black holes in AGN, GBHs and accreting 
neutron stars in their hard (low) states (see Tanaka & Lewin 1995; Gilfanov et al. 
1995 for the definition of the spectral states) show very similar X/7-ray spectra (see 
Zdziarski 1999; Barret et al. 2000). Furthermore, properties of their rapid temporal 
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variability are also similar (van der Klis 1995b; Wijnands & van der Klis 1999; 
Psaltis, Bclloni & van der Klis 1999: Ford et al. 1999; Edelson & Nandra 1999; 
Chiang et al. 2000). All this argues for a common mechanism of the X-ray production 
in all these sources. 

There are good reasons to believe that the main radiative mechanism for the 
production of the hard X-rays is Comptonization of soft photons (e.g., Shapiro, 
Lightman, & Eardley 1976; Sunyaev & Triimper 1979; Sunyaev & Titarchuk 1980). 
However, it is not completely clear what determines the observed spectral slopes. 
The geometry of the X-ray emitting region and the source of soft photons is still a 
matter of debate (see Svensson 1996; Poutanen 1998; Beloborodov 1999b; Wardzin- 
ski & Zdziarski 2000). 

An important clue to our understanding of the X-ray production came from the 
discovery of Fe lines (at <~ 6.4 kcV) and the hardening of the spectra above 10 keV 
in AGN (Pounds et al. 1990; Mushotzky, Done, & Pounds 1993; Nandra & Pounds 
1994), Cygnus X-l (e.g., Done et al. 1992; Gierlihski et al. 1997), and neutron stars 
(e.g., Yoshida et al. 1993). These features are associated with the reflection of hard 
X-rays from cold material (Basko, Sunyaev, & Titarchuk 1974; George & Fabian 
1991; Magdziarz & Zdziarski 1995; Poutanen, Nagendra, & Svensson 1996). These 
observations gave support to the so called two-phase accretion disc-corona models. 
In such models, X-rays are emitted by a hot rarified corona above the cold accretion 
disc (Haardt & Maraschi 1993; Haardt, Maraschi, & Ghisellini 1994; Stern et al. 
1995; Poutanen & Svensson 1996). Hard X-rays from the corona, being reprocessed 
in the cold disc, produce the reflection hump as well as most of the seed soft photons 
that are subsequently Comptonized to produce the hard X-rays. This is the feedback 
mechanism. The geometry of the corona determines the feedback factor which in its 
turn determines the spectral slope of the escaping radiation. The temperature of the 
emitting plasma (or to be more exact, the Kompanects y-parameter) is determined 
by the energy balance between heating (by magnetic reconnection?) and cooling 
(by Comptonization of soft photons). 

Further support for the feedback models was recently given by Zdziarski, Lu- 
binski, & Smith (1999) (see also Zdziarski 1999; Gilfanov, Churazov, & Rcvnivtsev 
2000) who found a correlation between the amount of reflection (R = f2/(27r), where 
Q is a solid angle subtended by cold material as viewed from the X-ray source) and 
the intrinsic photon spectral index, T, of the hard X-ray component. Such a correla- 
tion can easily be explained if there is overlap between the hot corona and the cold 
disc (Poutanen, Krolik, & Ryde 1997). The further the cold disc penetrates into 
the corona, the larger is the cooling, the smaller is the temperature of the corona, 
the softer is the spectrum, and, finally, the larger is the amplitude of the reflection. 
The model, however, appears to have trouble giving reflection amplitudes above 
-Rmax ~ 0.5 (if the coronal optical depth tt ~ 1, see Zdziarski et al. 1997) due to 
partial smearing of the reflection component by the hot corona. 

Alternatively, the observed R — T correlation can be reproduced by variations 
of the bulk velocity of the X/7-ray emitting plasma (Beloborodov 1999a, b). If the 
emitting regions arc sufficiently compact to produce electron-positron pairs, the 
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pressure of the radiation reflected and reprocessed in the disc accelerates pairs to 
mildly relativistic velocities away from the disc. On the other hand, for proton 
dominated plasmas, a small anisotropy in the energy dissipation mechanism can 
result in the ejection of particles away or towards the disc. Ejection away from the 
disc reduces R below 1 and leads to hard spectra, while ejection towards to the disc 
can result in apparent R > 1 as is observed in some objects. 

The physical possibility of the corona formation was studied by Galeev, Rosner, 
& Vaiana (1979). They showed that the magnetic fields, being amplified in the cold 
disc due to turbulent motions and differential rotation, do not have time to anni- 
hilate inside the disc on the inflow time scale. Instead, the field loops are expelled 
from the disc by buoyancy (the Parker instability) and they annihilate in the ten- 
uous corona. Beloborodov (1999a) showed that the mechanism studied by Galeev 
et al. is able to produce a corona of limited luminosity which is h/r (the ratio of 
the disc height to its radius) times smaller than the disc luminosity. By contrast, in 
some sources most of the energy escapes in the form of hard X-rays. Beloborodov 
also argued that the magneto-rotational instability (Velikhov 1959; Chandrasekhar 
1960; Balbus & Hawley 1991) increases the rate of the magnetic field generation 
(as compared with the Galeev et al. model) thus producing an active magnetic 
corona where a large fraction of the gravitational energy can finally be dissipated 
in magnetic flares. These qualitative arguments were recently supported by numeri- 
cal three-dimensional magnetohydrodynamical simulations of Miller & Stone (2000) 
who showed that about 25 % of the total energy dissipation can occur in the rarificd 
corona. 

An alternative to the magnetic corona is the hot disc model (Shapiro et al. 
1976; Ichimaru 1977; Narayan, Mahadevan, & Quataert 1998; Zdziarski 1998; Esin 
et al. 1998) which is also able to explain the observed X/7-ray spectra. In order to 
distinguish between the models, it would be helpful to compare the predictions of 
different models with the temporal variability data (see van der Klis 1995a,b and Cui 
1999a for recent reviews). Unfortunately, most of the papers on the spectral models 
do not consider the temporal variability. On the other hand, most of the models 
designed to explain the temporal variability data do not pay enough attention to 
the emission processes and the physics of the spectral formation. 

In this review, we will discuss the variability data keeping in mind recent ad- 
vances in modelling broad-band X/7-ray spectra of accreting black holes and neu- 
tron stars. Most attention will be paid to the time lags that can shed light on the 
mechanism of the X-ray production. Then we discuss simple phcnomcnological mod- 
els that are able to explain some of the observational facts. After that, we switch to 
the physical models. In particular, the properties of the Comptonizing regions will 
be discussed. We will point out the flaws in models that do not consider the energy 
balance in the "Compton cloud" , and then discuss models that satisfy the energy 
conservation law and confront them with the available data. 
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2. OBSERVING TIME LAGS IN ACCRETING BLACK HOLES AND NEU- 
TRON STARS 



The standard temporal characteristics that are usually computed are the power- 
density spectra (PDS) in different energy channels, auto/cross-correlation functions 
(ACF/CCF), the time/phase lags between the variability in different energy chan- 
nels, the coherence function, etc. Time lags have been studied by two methods, by 
constructing the CCF and by cross-spectral analysis (for details, see Lewin et al. 
1988; van der Klis 1989; Nowak et al. 1999a). 

2.1. Lags in Black Hole Sources 

The observations of Cygnus X-l from sounding rockets and HE AO 1 (Priedhorsky et 
al. 1979; Nolan et al. 1981) showed that the CCF between different energy channels 
peaks very close to zero lag (delay 40 ms), but it is slightly asymmetric. Similar 
asymmetry was found in the EXOSAT data by Page (1985) who claimed a 6± 1 ms 
shift of the peak of the CCF between the 5-14 keV and the 2-5 keV bands. Recent 
RXTE observations clearly show asymmetries of the CCFs, which, however, peak 
within ~ 1 ms from zero lag (see Fig. [I]). This suggests that the relation between the 
variation in the two bands are not simply a time delay. Asymmetry is also observed 
in the RXTE data of GX 339-4, where the CCFs are offset by < 5 ms from zero 
(using the 2-5 and 10-40 keV bands, Smith & Liang 1999). The CCFs of AGN also 
display similar properties (e.g., Papadakis & Lawrence 1995; Lee et al. 1999). 

The CCFs cannot be fitted with simple exponentials at any time scale. A 
reasonably good description of the CCF is in terms of a stretched exponential 
CCF{t) = /3exp[— (\t — tol/ry] (see the inset of Fig. |l|), where the normalisation 
w 1, the time where the CCFs peak to < 10~ 3 s, v ~ 2/3 in the range |i| < 0.3 s, 
and the time constant r is different for rising and decaying part of the CCF. Such a 
behaviour is probably the result of self-similarity of the light curve. It is interesting 
to note that the ACF of gamma-ray bursts also show a similar stretched exponential 
behaviour (see Stern & Svensson 1996; Beloborodov 1999c). 

Since the CCF does not show which frequencies contribute most to the observed 
lags, van der Klis et al. (1987) suggested to use instead the cross-spectrum for 
such an analysis.^] In the hard state of Cyg X-l observed by the Ginga satellite, 
the time lags between the variability in the 1.2-4.7 and 4.7-9.3 keV energy bands 
reached 0.1 s and had a strong Fourier-frequency dependence St(f) ~ i.e., 
the phase lag 8<f)(f) const (Miyamoto et al. 1988; Miyamoto & Kitamoto 1989). 
Similar lags were observed in other GBHs, GX 339-4 and GS 2023+338 (=V404 
Cyg), in their hard state (Miyamoto et al. 1992). The analysis of the RXTE data 
for Cyg X-l (Nowak et al. 1999a), GX 339-4 (Nowak et al. 1999b), IE 1740.7-2942 
and GRS 1758-258 (Smith et al. 1997), and GS 1354-644 (Revnivtsev et al. 2000) 

1 The cross-spectrum C(f) = S*(f)H(f), where S(f) and H(f) are the Fourier transforms of the 
light curves in the soft and hard energy channels, respectively. The phase lag, S<f>(f) = arg[C(/)], 
and the time lag, 5t(f) = 8tj>(f)/(2nf). The lags are positive when hard photons are lagging the 
soft ones. 
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FIGURE 1. The cross-correlation functions of Cyg X-l in the hard state (RXTE obser- 
vations from October 22, 1996). The solid curves are the autocorrelation function for the 
2-3.9 keV energy channel, the dotted and dashed curves are the CCFs for the 6-8.2 and 
14-70 keV vs. the 2-3.9 keV energy channel, respectively. Note, that CCFs in the hard state 
sometimes have much broader wings extending to ~ 8 s (e.g., Nolan et al. 1981; Maccarone 
et al. 2000). 

confirming the general features seen in the Ginga data, showed more complicated 
behaviour of the phase lag spectra which have a number of shelves and breaks (see 
Fig- 1). 

Grove et al. (1998) extended this analysis to higher energies with the data from 
CGRO/OSSFi. The time lags at low frequencies reached 0.3 s between the 50-70 
and the 70-100 keV photons in the light curve of GRO J0422+32 (=Nova Persei 

1992) . The breaks detected in the time lag spectrum at 0.1 Hz may be related to 
the quasi-periodic oscillation (QPO) observed at 0.23 Hz. The CGRO/BATSE data 
of Cyg X-l (Crary et al. 1998), GRO J0422+32 and GRO J1719-24 (=Nova Oph 

1993) (van der Hooft et al. 1999a, b) show very similar time lag spectra. 

The time lags of GBHs in their soft state (when T as 2.5) turn out to be somewhat 
different. Soft lags were observed between the 1.2-2.3 keV and the 2.3-4.6 keV bands 
in GX 339-4 (Miyamoto et al. 1991) and GS 1124-68 (=Nova Muscae 1991; see 
Miyamoto et al. 1993; Takizawa et al. 1997), while the higher energy photons were 
lagging the variability in the 2.3-4.6 keV band. In one observation of GS 1124-68, the 
variability in the 4.6-9.2 keV was the most advanced. The phase lag reached ~ 1 rad 
which is much larger than the lags in the broad band noise observed in the GBHs 
in their hard state. Rapid time variations were mostly due to the harder power-law 
component which is clearly seen from the rms amplitude. It is interesting that the 
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FIGURE 2. Time lags between signals in the 8.2-14 keV and the 3.9-6 keV bands vs 
the 2-3.9 keV band in Cyg X-l (RXTE observations from October 22, 1996). The dotted 
curves are the model of Poutanen & Fabian (1999b) with the maximum flare time scale of 
T m ax = 0.2 s. The solid curves show the same model with 70% of the observed Compton 
reflection produced at a distance of i re fl ~ 1 light seconds from the central X-ray source. 
The reflector acts as a low pass filter so that additional delays appear only at frequencies 
/ fs, 1/tteft, an d only at the energies where reflection is significant. See § 4.2 for details. 

lags seem to saturate above 10 keV. It is worth pointing out that the largest lags 
here are observed at the QPO frequency. On the other hand, the time lag spectrum 
of Cyg X-l in the soft state looks quite similar to that in the hard state (Cui et al. 
1997). 

Time lags have been observed in GRS 1915+105 in the 67 mHz QPO by Cui 
(1999b) and in the broad-band noise and QPOs by Reig et al. (2000). The lags 
show a very complicated structure, sometimes changing signs from one harmonic 
to another, and the sign also depends on the frequency of the QPO. Wijnands, 
Homan, & van der Klis (1999) and Cui, Zhang, & Chen (2000) observed lags in the 
broad-band noise and the QPOs of XTE J1550-564. 

2.2. Lags in Neutron Star Sources 

Hard time lags in neutron star sources were discovered by Hasinger (1987) in the 
CCF of Cyg X-2 (comparing the 1-5 and the 5-17 keV bands) in its horizontal branch 
(HB, see Lewin et al. 1988 for definitions of branches). The CCF also had sinusoidal 
oscillations due to the QPO. The time lags showed anti-correlation with the QPO 
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frequency, dropping from 4 ms to 1.5 ms when the QPO changed from 20 Hz to 50 
Hz. Hasinger interpreted the lags as delays due to scattering (Comptonization) in 
the hot cloud and the anti-correlation as an indication of a change in the system 
size. Associating the QPO frequency with the Keplerian frequency at some radius 
gives the relation, St(f) oc Jq P q, while the actual data are much better described 
by St oc fqp , i.e. 8(f) = 2-rrf ■ St — const. 

Using the cross-spectrum techniques van der Klis et al. (1987) confirmed the 
existence of - 3 ms hard lags in the 20-40 Hz QPO of Cyg X-2 (and GX 5-1) 
and discovered 8 ms soft lags in the low-frequency noise, which were interpreted as 
a softening of the spectrum during the shots that cause QPOs. These results were 
confirmed by Vaughan et al. (1994) who also showed (from the analysis of the Ginga 
data of GX 5-1 on the HB) that the time lags increase with photon energy. 

In both Cyg X-2 and GX 5-1 in their normal branch, the lags at the ~ 5 Hz 
QPO showed energy dependence (Mitsuda & Dotani 1989; Vaughan et al. 1999) 
reaching 5<p ~ 7r rad (i.e., St ~ 0.2s) for 10 keV photons vs 2 keV. At the same 
time, the rms amplitude of the QPO in Cyg X-2 had a minimum at 5 keV and in 
GX 5-1 it increased above 2.5 keV. This behaviour can be interpreted as a pivoting 
of the spectrum around 3-5 keV. 

With the larger effective area of RXTE, Ford et al. (1999) and Olive & Bar- 
ret (2000) discovered phase lags in the broad-band noise of three atoll sources, 
4U0614+09, 4U1705-44, and 4U1728-34. These lags are very similar to those in 
GBHs like Cyg X-l and GX 339-4, which tells us that the mechanism responsible 
for the lags does not depend on the presence or absence of the hard surface of the 
neutron star, magnetosphere, boundary layer, etc., but instead is a property of the 
accretion flow. 

A number of neutron stars show kHz QPOs in their light curves as revealed 
by RXTE. Kaaret et al. (1999) find 25 fis soft lags between the 4-6 keV and the 
> 9 keV photons in the 800 Hz QPO in the atoll source 4U1636-536. Analysing 
the 550 Hz oscillations of Aquila X-l, Ford et al. (1999) found soft lags, 8<j) ~ 1 
rad, between the 3-6 keV and the > 6 keV photons. Similar lags were found in the 
accreting millisecond pulsar SAX J1808.4-3658 (Cui, Morgan, & Titarchuk 1998; 
Ford 2000). The lags in the 830 Hz QPO in 4U 1608-52 (Vaughan et al. 1997, 1998) 
reach 60 /zs between the 5 and the 25 keV photons. 

3. PHENOMENOLOGICAL ENERGY DEPENDENT SHOT NOISE MODELS 

It is clear that the observed zoo of time lags cannot be explained by any single 
model. Different mechanisms should be involved in producing the lags at different 
Fourier frequencies, the lags in the QPOs and the coherent pulsations, and the lags 
in the broad-band noise. 

Let us first consider the simplest possible model that produces time lags: a shot 
noise model (Terrell 1972), where shots (=flares) are uncorrelated with each other. 
We assume that the shot time profiles at different energies have the same shape, 
but slightly different time constants. As an example we take a shot profile at soft 
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energies s(t) = [t/r] p cxp[— t/r] and at hard energies h{t) = [t/(r]T)) p exp[— t/(t]T)], 
where t > is measured from the beginning of the shot and p is positive. The 
Fourier transforms, S(f) and H(f), are S(f) cx tT( P + 1)/(1 - i27r/r) p+1 and 
H(f) cx r]TT{p+l)/(l~i2TTfriT)P +1 . ThePDSsarecx \S{f)\ 2 and \H(f)\ 2 . For small 
frequencies, / <§; 1/(2itt), the PDSs have a flat dependence on frequency, cx f°, 
while for large frequencies, / ^S> l/(2nr), the PDSs decay as f~ 2 (j>+ 1 ). The power 
per logarithm of frequency (i.e., / x PDS(f)) peaks for soft photons at f s ,max = 
l/[2irTy/2p + 1J and at fh,max = l/[r]2nTy / 2p~+T\ for hard photons. The phase lags 
8<j)(f) = (p + l)[arctan(?727r/r) — arctan(27r/r)]. The lag is positive when hard 
photons are lagging soft ones (i.e., for 77 > 1). For small frequencies, Scf>(f) rises as 
ps (p+l)(ry— 1)27t/t, while for large frequencies, it decays as (j>+l)(r}—l)/(r]2irfT). 
The lag reaches a maximum of 8<f> max = 2(p+l)(arctany / ?7— 7r/4) at / = \j(2-KT^ff\) 
close to /.j, max and max , the frequencies where / x PDS s ,h(f) peak. 

One can also consider a modified shot noise model, where the shot time scales are 
distributed according to a power law, p(r) cx t~ p between r m ; n and r max (see, e.g., 
Miyamoto & Kitamoto 1989; Lochner, Swank, & Szymkowiak 1991), with the same 
ratio r). Physically this could correspond to, for example, the situation when flares 
of different durations appear at different radii from the central object (Poutanen & 
Fabian 1999a). A power-law distribution of r assures that the PDS is also a power- 
law cx / _ ( 3_ p) (Lochner et al. 1991). If the flares are self-similar, then the phase lag 
will be constant ps <5</> max for / > / min = l/(27rr max ) and / < / max = l/(27rr min ), 
decay as 1// at / > / max and rise linearly at / < / m i n - The corresponding time lags 
are constant, dt max = 27TT max <5</> max , for / < / m i n , and decay approximately as 1// 
between / m i n and / max - Note that in this model the coherence function (Vaughan & 
Nowak 1997; Nowak et al. 1999a) is close to unity, since the light curves at different 
energies are almost perfectly synchronised. 

If we assume that r/ > 1, there are hard lags and the predicted behaviour of the 
time lags and coherence function is in a very good agreement with the observations 
of GBHs (Poutanen & Fabian 1999a). However, this model contradicts the CCF of 
Cyg X-l (see Fig. | and Maccarone et al. 2000). The CCF becomes narrower at 
larger energies which requires the shots to be narrower at larger energies. If one, 
however, reverses the time profiles of the shots, so that they rise slower and decay 
faster (e.g., s(t) = (— t/r) p exp(i/V), t < 0), and one assumes ij < 1 (i.e., hard shots 
are narrower), the CCFs and the time lags can be reproduced, simultaneously. 

Much more complicated models which also account for lags in the QPOs sources 
were developed by Shibazaki et al. (1988). We just note here that if a signal consists 
of shots appearing almost periodically and if shots at different energies are shifted 
in time one against another (or, e.g., the minima are reached at the same time and 
the peaks are not), the phase lag has a very complicated dependency on frequency 
(e.g., changes sign from one harmonic to another) depending on the shot profiles. 
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4. PHYSICAL MECHANISMS FOR PRODUCING LAGS 



4.1. Static Compton Cloud Models 

Since Comptonization is the most probable mechanism for X-ray production in 
compact objects, it is natural to attribute the time delays between hard and soft 
photons to this process. Hard photons arc the result of more scattering and so 
emerge after, or lag behind, softer ones. Consider a static "Compton cloud" with 
fixed Thomson optical depth, tt, and electron temperature O = kT e /m e c 2 . A soft 
seed photon of energy E injected into the cloud increases its energy by a factor of 
Ai = 1 + 40 + 160 2 on average after each scattering, so that after TV-scatterings 
its energy = E . The photon mean free path is A w Rj max(l, tt) (where R 
is the size of the X-ray producing region, and where we accounted for the fact that 
we are interested only in those photons that actually have undergone scatterings in 
the cloud). The time between successive scatterings is then t c = i?/(c max[l, tt]), 
so the time needed to reach the energy Em is (Sunyaev & Titarchuk 1980; Payne 
1980) 

au R l c HEn/E ) 

t N = iVi c = — - — , 

max(l,Tx) m^li 

which translates to t N ~ 10~ 4 s for kT e ~ 50 keV, tt ~ 1, R = 10 km, and 
E N /E ~ 10. 

This model was criticised by Miyamoto et al. (1988), Miyamoto ct al. (1991) 
and Vaughan et al. (1994). First, the large size of the cloud (10 3 — 10 5 i? g , where R g 
is the Schwarzschild radius, 2GM/c 2 ) is needed to produce large delays observed 
in GBHs and neutron stars. Such cloud is physically unrealistic, since most of the 
gravitational energy is dissipated within 10R g . Second, the lags predicted by the 
model are independent of the Fourier frequency (Miyamoto et al. 1988) contrary 
to the observed ~ 1/f dependence. Third, it is assumed that the soft photons 
produce the variability, while the hot cloud is not variable. Observationally, it is 
well established that when a soft black body spectrum is observed in GBHs it is 
much less variable than the hard X-rays (e.g. Miyamoto et al. 1991), so that the 
hard X-ray variability is most probably intrinsic to the hot cloud itself. 

Finally, we would like to point out that due to the requirement of energy con- 
servation the whole concept of a static Compton cloud with a constant temperature 
is physically unrealistic. The total emitted X-ray luminosity (produced by Comp- 
tonization of soft radiation) is L to t = Lh + L Sl where Lh is the heating rate in 
the hot cloud and L s is the luminosity of seed soft photons. For hard spectra (i.e. 

Ltot ~ Lh- The total X-ray luminosity is thus a function of the heating 
rate only and it does not depend on the amount of seed soft photons. By chang- 
ing L s , one effectively changes the spectral slope of the emergent X-ray radiation 
which is a function of the Compton amplification factor A = Lh/L s . This results 
in the pivoting of the spectrum (see Poutanen 1998; Beloborodov 1999b) without a 
noticeable increase in Ltot- The larger the L s , the smaller the equilibrium temper- 
ature of the emitting electrons, and the softer the spectrum. By contrast, in static 
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Compton cloud models no changes in the electron temperature are considered in 
reaction to the changes in the number of soft photons, violating thus the energy 
conservation law. In order to increase the emitted luminosity, one has to change the 
energy dissipation rate in the cloud, but then exactly these changes will be driving 
the variability. 

Recently, Kazanas, Hua, & Titarchuk (1997) (see also Bottcher & Liang 1998; 
Hua, Kazanas, & Cui 1999) modified the Comptonization model. Instead of a ho- 
mogeneous Compton cloud, the density profile, n(r) ~ 1/r, was assumed (in this 
case, one has equal optical depth per logarithm of radius). The variability is still 
driven by changing the rate of soft photon injection in the center of the cloud. Then 
larger radii produce lower frequency variability (filtering out high frequency signal) 
and larger lags, while the smaller radii produce higher frequency variability and 
smaller lags (see also Nowak et al. 1999a, c). This model solves only one of the prob- 
lem mentioned above (1// time lag dependence), while the other problems remain 
unsolved. Another modification of the model was considered by Bottcher & Liang 
(1999) based on an earlier suggestion by Miyamoto et al. (1988). Here, small cold 
clouds are assumed to free-fall into the hot central cloud thus changing the input 
of soft photons. Again, this model has problems with energy conservation. 

4.2. The Dynamic Compton Cloud 

Miyamoto et al. (1988) pointed out that some modulation mechanism must be 
invoked to produce the strongly frequency-dependent time lags. Poutanen & Fabian 
(1999a, b) proposed a model where the time lags are produced by the evolution of 
the flare spectrum. They assumed that the energy dissipation varies in time. For 
a small (of the order of R g ) emitting region (ER) one can consider the spectral 
evolution as a sequence of steady-states as long as the characteristic time scale 
of variability is r ^ R/c. Any changes in the amplification factor A would cause 
spectral variability and, specifically, a continuous increase of A with time during 
the course of the flare would cause a soft-to-hard spectral evolution producing hard 
time lags of the order of the flare time scale (see § ^]). Poutanen & Fabian (1999b) 
considered three mechanisms that can increase A. 

(1) A flare starts in the background of soft photons. The X-ray spectrum is soft 
as long as Lh < L s ,bkg- With increasing Lh, the soft photon input gets dominated 
by reprocessed photons. The spectral slope is then determined by the feedback 
parameter D i=s l/A (the fraction of Lh returned to the ER after the reprocessing in 
the disc into soft seed photons, see Stern et al. 1995; Svensson 1996; Beloborodov 
1999b) and the spectrum becomes hard. 

(2) The dissipation is accompanied by pumping net momentum into the hot plasma 
of the emission region (Beloborodov 1999a,b). The resulting bulk velocity increases 
with increasing luminosity. It leads to a lower feedback and higher A (if the velocity 
is directed away from the disc). 

(3) The differential rotation of the footpoints of a magnetic loop at the disc surface 
causes a twisting and elevation of the loop (Romanova et al. 1998). The time scale of 
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the evolution is of the order of the Keplerian time-scale. When the emission region 
moves away from the disc, the feedback decreases and A increases (see Fig. ||). 

In all these cases the spectral evolution proceeds from soft to hard. The hard 
time lags between energies E and Eq are oc t\h{E/Eq). If there is a distribution 
of time scales r between, say, 1 ms and 0.3 s (e.g., Keplerian time scales at the 
radii between the innermost radius of the accretion disc and ~ 50R g ), the time lags 
oc 1// at the characteristic frequencies of the variability (see § [}] and Fig. ||). 

If t ~ a few light crossing time of the ER, the spectrum in the beginning of the 
flare is hard because of photon starvation (one needs a few R/c to get reprocessed 
soft photons into the ER) and softens towards the end of the flare (Poutanen & 
Fabian 1999a; Malzac & Jourdain 2000) producing soft time lags. Observing the 
change of sign of the time lags at some Fourier frequency, / sg n, would determine 
the size of the ER that produces variability at these frequencies (e.g., for Cyg X-l, 
R< 10i? 9 (30 Hz// Sgn )). 




E (kcV) Time, t/r 

FIGURE 3. (a) Spectral evolution of a magnetic flare. Time resolved spectra (without 
Compton reflection) are presented by dotted curves; marks are times, t/r, from the be- 
ginning of the flare. The time-averaged Comptonized spectrum is shown by a thin solid 
curve. The time-averaged spectrum of Cygnus X-l (simultaneous Gvnga and OSSE data 
from June 1991, the data set # 1 in Gierlihski et al. 1997) is plotted with crosses and the 
best fit with the flare model with a solid curve (x 2 /dof= 50.0/75). Interstellar absorption 
is removed when plotting the model spectrum, (b) The flare light curves at 3, 9, 27, and 81 
keV are presented by dotted curves. Solid curve - the heating rate, lh(t) oc (t/r) 2 exp(— t/r); 
dashed curve - the feedback factor D(t); dot-dashed curve - the temperature of the emission 
region. See Poutanen & Fabian (1999a,b) for details. 
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4.3. Small Scale Spectral Transitions 

The models considered above can explain the time lags in the broad-band noise. 
In some QPO sources the time lags show a very complicated behaviour (see § ^) 
and may require different explanation. In radiation- hydrodynamic model (see, e.g., 
Lamb 1989; Miller & Lamb 1992), QPOs appear as a result of oscillations in the 
optical depth of the radial flow due to the radiation feedback from the neutron star 
surface. The resulting spectral pivoting produces phase lags and the increase of the 
rms amplitude variability above the pivoting point. This model, however, is not 
applicable to the QPOs and the lags in black holes sources (see, e.g., Takizawa et 
al. 1997) because of the absence of a hard surface. 

The galactic microquasar GRS 1915+105 shows large amplitude oscillations with 
periods varying from <1 up to 100 s. These time scales are a few orders of mag- 
nitudes larger than Keplerian time scales and up to 10 6 time larger than the light 
crossing time of the ER. The best candidate for producing the spectral variability 
that causes the lags in GRS 1915+105 is the oscillation of the inner radius of the 
cold disc on viscous time scales. Such oscillations are similar to the spectral tran- 
sitions, but have smaller amplitude and occur at shorter time scales. Changes of 
the relative geometry of the hot corona and the cold disc (with or without changes 
of the total luminosity) cause spectral pivoting at a few keV (see, e.g., Poutanen, 
Krolik, & Ryde 1997; Esin et al. 1998). The fluxes below and above the pivot point 
oscillate then with a phase shift of ~ tt. The rms amplitude of the QPO increases 
with the energy. The phase lags between the energies above the pivot point can 
then be produced if the oscillations are time asymmetric (see Fig. 11 in Morgan, 
Remillard, & Greiner 1997; Vilhu & Nevalainen 1998). 

Similar (but aperiodic) changes in the inner disc radius can be responsible for the 
broad-band variability observed at / < 1 Hz in, e.g., Cyg X-l. Associated spectral 
changes can manifest themselves in time lags observed at these frequencies. 



4.4. Delays due to Compton Reflection and Reprocessing 

The spectra of accreting GBHs and neutron stars show signatures of Compton 
reflection (see § [[]). Some fraction of the X-ray photons can be reflected from the 
outer edge of a flared accretion disc, a wind from the companion, etc. Such a reflector 
acts as a low pass filter smearing out the high frequency variations and produces 
lags corresponding to the light travel time to the reflector only at lower frequencies. 
Such processes can explain the break in the time lag spectra observed in Cyg X-l 
(see Fig. §) and in other GBHs at / < 1 Hz. 

Reprocessed soft radiation which accompanies Compton reflection is emitted 
in the optical and UV spectral bands if the reprocessing occurs far away from the 
central X-ray source. The time delays can then be measured between the optical/UV 
and the X-ray radiation (e.g., Hynes et al. 1998). On the other hand, reprocessing 
in the vicinity of the X-ray emitting region, produces time delays of the order of the 
time scale of the spectral evolution in the hard X-ray band. This may be one reason 



for the observed soft lags in the soft state of GX 339-4 and GS 1124-68 (§ |2.1.|) 
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4.5. Hot Spots on the Neutron Star Surface 

Some neutron star sources show lags in their periodic oscillations at kHz frequencies. 
Ford et al. (1999) and Ford (2000) interpreted the soft lags in Aquila X-l and in the 
accreting millisecond pulsar SAX 1808.4-3658 using a model of a rotating hot spot 
with a black body spectrum at the surface of a neutron star where the lags appear 
due to Doppler effects. The weak energy dependence of the rms amplitude of the 
oscillations reported by Cui et al. (1998) rules out the black body model for the 
spectrum used by Ford. Detailed analysis of the pulsations in the time domain by 
folding techniques by Revnivtsev (1999) revealed that the pulse profile is distorted 
at different energies, while the minima are reached at the same time (i.e., there are 
no lags in the normal meaning of this word) . 

5. SUMMARY 

Time lags and other temporal variability data provide strong constraints on the 
models of the X-ray production. It was demonstrated that static Compton cloud 
models are based on physically unrealistic assumptions. The models invoking spec- 
tral evolution of the flare spectrum can fit both the CCF and the time lag Fourier 
spectra only if (1) the energy dissipation rate increases slowly and decreases rapidly 
and (2) the flare spectrum evolves from soft to hard. If soft seed photons are pro- 
duced by reprocessing the hard ones, the change of sign in the time lag spectrum is 
expected at high frequencies corresponding to the light crossing time of the emis- 
sion region. The absence of such a change would put constraints on the size of the 
emitting region. 

We also argued that the reflection of hard X-rays from the outer part of the 
accretion disc produces time delays that we already might have observed in GBHs. 
If so, the disc should be flared and the break in the time lag Fourier spectra then 
corresponds to the size of the accretion disc. Of course, such an interpretation is 
not unique. Alternatively, small scale spectral transitions (e.g., oscillations of the 
inner radius of the accretion disc at viscous time scales) might produce time lags 
observed at lower frequencies. 

In the case of (quasi-) periodic oscillations from the neutron star sources, we 
argued that in order to reproduce both the time lags and the energy dependent rms 
amplitude, the spectrum of the hot spots should not be close to a black-body. 

ACKNOWLEDGEMENTS 

This work was supported by the Swedish Natural Science Research Council and the 
Anna-Greta and Holgcr Crafoord Fund. I thank Katja Pottschmidt for providing 
the time lag Fourier spectra and the light curves of Cyg X-l used in the calculations 
of the cross-correlation functions. I am grateful to Andrei Beloborodov and Roland 
Svensson for valuable comments. 



13 



REFERENCES 

Balbus, S.A., Hawley, J.F. 1991, ApJ, 376, 214 

Basko, M.M., Sunyaev, R.A., Titarchuk, L.G. 1974, A&A, 31, 249 

Barret, D., Olive, J.F., Boirin, L., Done, C, Skinner, G.K., Grindlay, J.E. 2000, ApJ, in press 
jastro-ph/991104S ) 



Bcloborodov, A.M. 1999a, ApJ, 510, L123 

Bcloborodov, A.M. 1999b, in High Energy Processes in Accreting Black Holes, ASP Conf. Series 



Vol. 161, ASP, San Francisco, p. 295 (astro-ph/990110j ) 



Beloborodov, A.M. 1999c, in Gamma-Ray Bursts: The First Three Minutes, ASP Conf. Series Vol. 



190, ASP, San Francisco, p.47 ( istro-ph/9911122) 
Bottcher, M., Liang, E.P 1998, ApJ, 506, 281 
Bottcher, M., Liang, E.P. 1999, ApJ, 511, L37 
Chandrasekhar, S. 1960, Proc. Natl. Acad. Sci. USA, 46, 253 
Chiang, J. et al. 2000, ApJ, 528, 292 
Crary, D.J. et al. 1998, ApJ, 493, L71 

Cui, W., Zhang, S.N., Focke, W., Swank, J.H. 1997, ApJ, 484, 383 
Cui, W., Morgan, E.H., Titarchuk, L.G. 1998, ApJ, 504, L27 

Cui, W. 1999a, in High Energy Processes in Accreting Black Holes, ASP Conf. Series Vol. 161, 



ASP, San Francisco, 97 ( ^stro-ph/9809408| ) 
Cui, W. 1999b, ApJ, 524, L59 

Cui, W., Zhang, S.N., Chen, W. 2000, ApJ, 531, L45 

Done, C, Mulchaey, J.S., Mushotzky, R.F., Arnaud, K.A. 1992, ApJ, 395, 275 
Edelson, R., Nandra, K. 1999, ApJ, 514, 682 

Esin, A.A., Narayan, R., Cui, W., Grove, E.C., Zhang, S.-N. 1998, ApJ, 505, 854 
Ford, E.C. et al. 1999, ApJ, 512, L31 



Ford, E.C. 2000, ApJ, submitted ( |astro-ph/0002052 ) 
Galeev, A.A., Rosner, R., Vaiana, G.S. 1979, ApJ, 229, 318 
George, I.M., Fabian, A.C. 1991, MNRAS, 249, 352 
Gierlihski, M. et al. 1997, MNRAS, 288, 958 

Gilfanov, M. et al. 1995, in The Lives of the Neutron Stars, NATO C 450. Kluwer Academic 
Publishers, Dordrecht, p. 331 



Gilfanov, M., Churazov, E., Revnivtsev, M. 2000, A&A, in press (astro-ph/9910084) 
Grebenev, S.A. et al. 1993, ApJS, 97, 281 

Grebenev, S.A., Sunyaev, R.A., Pavlinsky, M.N. 1997, Adv. Space Res., 19, (1)15 

Grove, J.E. et al. 1997, in Proceedings of 4th Compton Symposium, AIP Conf. Proc. Vol. 410, 

AIP, New York, p. 122 
Grove, J.E. et al. 1998, ApJ, 502, L45 
Haardt, F., Maraschi, L. 1993, ApJ, 413, 507 
Haardt, F., Maraschi, L., Ghisellini, G. 1994, ApJ, 432, L95 

Hasinger, G. 1987, in The Origin and Evolution of Neutron Stars, IAU Symp. 125, D. Reidel Publ. 

Co., Dordrecht, 333 
Hua, X.-M., Kazanas, D., Cui, W. 1999, ApJ, 512, 793 

Hynes R.I., O'Brien, K., Home, K., Chen, W., Haswell, C.A. 1998, MNRAS, 299, L37 



14 



Ichimaru, S. 1977, ApJ, 214, 840 

Johnson, W.N. et al. 1997, in Proceedings of 4th Compton Symposium, AIP Conf.Proc. Vol. 410, 

AIP, New York, p. 283 
Kaaret, P., Piraino, S., Ford, E.C., Santangelo, A. 1999, ApJ, 514, L31 
Kazanas, D., Hua, X.-M., Titarchuk, L. 1997, ApJ, 480, 735 

Lamb, F.K. 1989, in Proc. 23rd ESLAB Symp. on Two- Topics in X-ray Astronomy, ESA SP-296, 
p.215 

Lee, J.C. et al. 1999, MNRAS, submitted ( astro-ph/9909239 ) 

Lewin, W.H.G., van Paradijs, J., van der Klis, M. 1988, Space Sci. Rev., 46, 273 

Lochner, J.C, Swank, J.H., Szymkowiak, A.E. 1991, ApJ, 376, 295 

Maccarone, T. et al. 2000, in preparation 

Magdziarz, P., Zdziarski, A. A. 1995, MNRAS, 273, 837 

Malzac, J., Jourdain, E. 2000, A&A, submitted 

Miller, G.S., Lamb, F.K. 1992, ApJ, 388, 541 

Miller, K.A., Stone, J.M. 2000, ApJ, in press ( astro-ph/9912135 ) 

Mitsuda, K., Dotani, T. 1989, PASJ, 41, 557 

Miyamoto, S., Kitamoto, S. 1989, Nature, 342, 773 

Miyamoto, S., Iga, S., Kitamoto, S., Kamado, Y. 1993, ApJ, 403, L39 

Miyamoto, S., Kitamoto, S., Mitsuda, K., Dotani, T. 1988, Nature, 336, 450 

Miyamoto, S. et al. 1991, ApJ, 383, 784 

Miyamoto, S. et al. 1992, ApJ, 391, L21 

Morgan, E.H., Remillard, R.A., Greiner, J. 1997, ApJ, 482, 993 

Mushotzky, R.F., Done, C., Pounds, K.A. 1993, Ann. Rev. Astron. Astrophys., 31, 717 
Nandra, K., Pounds, K.A. 1994, MNRAS, 268, 405 

Narayan, R., Mahadevan, R., Quataert, E. 1998, in Theory of Black Hole Accretion Discs, Cam- 
bridge Univ. Press, Cambridge, p. 148 
Nolan, P.L. et al. 1981, ApJ, 246, 494 
Nowak, M.A. et al. 1999a, ApJ, 510, 874 
Nowak, M.A., Wilms, J., Dove, J.B. 1999b, ApJ, 517, 355 
Nowak, M.A. et al. 1999c, ApJ, 515, 726 
Olivie, J.-F., Barret, D. 2000, these proceedings 
Page, C.G. 1985, Space Sci. Rev., 40, 387 
Papadakis, I.E., Lawrence, A. 1995, MNRAS, 272, 161 
Payne, D.G. 1980, ApJ, 237, 951 
Pounds, K.A. et al. 1990, Nature, 344, 132 

Poutancn, J. 1998, in Theory of Black Hole Accretion Discs, Cambridge Univ. Press, Cambridge, 
p.100 

Poutanen, J., Fabian, A.C. 1999a, MNRAS, 306, L31 

Poutanen, J., Fabian, A.C. 1999b, in High Energy Processes in Accreting Black Holes, ASP Conf. 

Series Vol. 161, ASP, San Francisco, p. 135 
Poutanen, J., Svensson, R. 1996, ApJ, 470, 249 
Poutanen, J., Krolik, J.H., Ryde, F. 1997, MNRAS, 292, L21 
Poutanen, J., Nagendra, K.N., Svensson, R. 1996, MNRAS, 283, 892 
Priedhorsky, W. et al. 1979, ApJ, 233, 350 



15 



Psaltis, D., Belloni, T., van der Klis, M. 1999, ApJ, 520 , 262 
Reig, P. et al. 2000, ApJ, submitted (astro-ph/0001134 1 

Revnivtsev, M. 1999, PhD thesis, Space Research Institute, Moscow (astro-ph/9912556) 
Revnivtsev, M., Borozdin, K., Priedhorsky, W.C., Vilhlinin, A. 2000, ApJ, 530, in press (astro 



ph/9905380) 



Romanova, M. M. et al. 1998, ApJ, 500, 703 

Shapiro, S.L., Lightman, A.P., Eardley, D.N. 1976, ApJ, 204, 187 
Shibazaki, N. et al. 1988, ApJ, 331, 247 
Smith, D.M. et al. 1997, ApJ, 489, L51 
Smith, I. A., Liang, E.P. 1999, ApJ, 519, 771 

Stern, B.E., Poutanen, J., Svensson, R., Sikora, M., Begelman, M.C. 1995, ApJ, 449, L13 

Stern, B.E., Svensson, R. 1996, ApJ, 469, L109 

Sunyaev, R.A., Titarchuk, L.G. 1980, A&A, 86, 121 

Sunyaev, R.A., Trumper, J. 1979, Nature, 279, 506 

Svensson, R. 1996, A&AS, 120C, 475 

Takizawa, M. et al. 1997, ApJ, 489, 272 

Tanaka, Y., Lewin, W.H.G. 1995, in X-ray binaries, Cambridge Astrophysics Series, vol. 26, Cam- 
bridge University Press, Cambridge, 126 
Terrell, N. J. Jr. 1972, ApJ, 174, L35 
van der Hooft, F. et al. 1999a, ApJ, 513, 477 
van der Hooft, F. et al. 1999b, ApJ, 519, 332 
van der Klis, M. et al. 1987, ApJ, 319, L13 

van der Klis, M. 1989, in Timing Neutron Stars, NATO ASI C 262. Kluwer Academic Publishers, 
27 

van der Klis, M. 1995a, in The Lives of the Neutron Stars, NATO ASI C 450. Kluwer Academic 
Publishers, 301 

van der Klis, M. 1995b, in X-ray binaries, Cambridge Astrophysics Series, vol. 26, Cambridge 

University Press, Cambridge, p. 252 
Vaughan, B.A., Nowak, M.A. 1997, ApJ, 474, L43 
Vaughan, B. et al. 1994, ApJ, 421, 738 

Vaughan, B. et al. 1997, ApJ, 483, L115 (erratum 1998, ApJ, 509, L145) 
Vaughan, B.A. et al. 1999, A&A, 343, 197 
Velikhov, E.P. 1959, Sov. Phys. JETP, 36, 995 
Vilhu, O., Nevalainen, J. 1998, ApJ, 508, L85 

Wardzinski, C, Zdziarski, A. A. 2000, MNRAS, in press ( astro-ph/9911126 ) 

Wijnands, R., van der Klis, M. 1999, ApJ, 514, 939 

Wijnands, R., Homan, E., van der Klis, M. 1999, ApJ, 526, L33 

Yoshida, K. et al. 1993, PASJ, 45, 605 

Zdziarski, A. A. 1998, MNRAS, 296, L51 

Zdziarski, A. A., Johnson, W.N., Poutanen, J., Magdziarz, P., Gierlihski, M. 1997, in The Trans- 
parent Universe, ESA SP-382, p.373 

Zdziarski, A. A. 1999, in High Energy Processes in Accreting Black Holes, ASP Conf. Series Vol. 
161, p.16 

Zdziarski, A.A., Lubinski, P., Smith, D.A. 1999, MNRAS, 303, Lll 



1G 



